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Abstract. - 

The problem of the occurrence of an energy cascade for Alfvenic turbulence in solar wind 
plasmas was hystorically addressed by using phenomenological arguments based to the weak- 
ness of nonlinear interactions and the anisotopy of the cascade in wave vectors space. Here, this 
paradox is reviewed through the formal derivation of a Yaglom relation from anisotropic Mag- 
netohydrodynamic equation. The Yaglom relation involves a third-order moment calculated 
from velocity and magnetic fields and involving both Elsasser vector fields, and is particularly 
useful to be used as far as spacecraft obervations of turbulence are concerned. 



Introduction. - In a seminal paper, Dobrowolny, Mangeney and Veltri ^^ (hereafter 
DMV), rised the question of the existence of an energy cascade in Magnetohydrodynamic 
(MHD) turbulence because the apparent contradiction of two competing observations within 
the solar wind turbulence by in situ satellite measurements. In fact, since the oldest space 
flights, both a well defined turbulent spectrum, and strong correlations between velocity and 
magnetic field fluctuations have been observed [2113^ (for a modern review of turbulence in 
the solar wind cfr. Ref. [3]). The apparent contradiction between these observations can be 
immediately seen by introducing the Elsasser variables 



where and represent the velocity and magnetic field respectively, while p is the mass 
density. These quantities represent Alfvenic fiuctuations propagating along the background 
magnetic field, in opposite directions. MHD equations can be immediately written in terms 
of these variables as 



where tt = P/ p (P being the the total pressure), dt represents time derivative while di rep- 
resents derivative with respect to the spatial variable Xi. The kinematic viscosity v and the 
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magnetic diffusivity n form the dissipative coefficients = (t^ ± /i)/2. The second term in 
equations ([1]) shows that nonhnear interactions only occur between opposite sign fluctuations. 
Since high correlations between velocity and magnetic fluctuations imply either = or 
Zf^ = 0, a turbulent energy cascade should be incompatible with the disappearence of one of 
the alfvenic fluctuations. The puzzle have been apparently solved by DMV Ij. In presence of 
a strong magnetic field, nonlinear interactions are slowed down by the transport of fluctua- 
tions (Alfven effect). The usual Kolmogorov's phenomenlogy must then be modified in favor 
of the Iroshnikov-Kraichnan (IK). This yields to the fact that the energy transfer rates per 
unit mass for both pseudo-energies associated to alfvenic fluctuations must be of the same 
order, e"*" ~ e~. More precisely, they must have the same scaling laws in the IK phenomenol- 
ogy [7]. Thus, an initial small unbalance between alfvenic fluctuations is maintained during 
the cascade, eventually leading to both a turbulent spectrum, and high correlations P^. This 
framework is referred to in the literature as Alfvenic turbulence. The above arguments have 
been criticized [51[3] on the basis of the fact that, at variance with the conjecture in Ref. [I], 
both in closures equations [5| and in direct numerical simulations [51 the energy transfer rates 
are never the same. A different phcnomenological argument has also been introduced 0. In a 
strongly anisotropic medium as the solar wind, the energy cascade will eventually develop in 
the direction transverse to the background magnetic field. Using the hypothesis of a critical 
balance between the time of transport in the parallel direction and the eddy turnover time for 
turbulence in the transverse direction, a Kolmogorov's spectrum, rather than a Kraichnan's 
spectrum, is expected for transverse wavevectors. 

In homogeneous and isotropic fluid turbulence, the energy cascade is evidenced by the 
observation of a well defined relation between the third-order longitudinal structure function 
and the energy dissipation transfer rate, namely the well known 4/5-Kolmogorov's law ■9j. 
A similar relation have been derived for MHD, following the Yaglom law for passive scalars, 
in the framework of homogeneous and isotropic turbulence jlOj . This law has been recently 
observed in samples of ecliptic [12J and polar jll) solar wind, thus showing unambiguously 
that a turbulent cascade is active. These observations solve the apparent paradox rised in 
the DMV paper, confirming the presence of both strong correlations and turbulent cascade. 
However, the anisotropic nature of solar wind turbulence occasionally violates the conditions 
for the observation of the Yaglom law, namely of the cascade. In this letter we analyze the 
relevance of anisotropy conjecture in deriving the Yaglom's scaling law for MHD turbulence. 

The Yaglom law for anisotropic and isotropic MHD. - Consider the anisotropic MHD 
equations ([T} written twice for Elsasser variable Z^{xi) at the point Xi, and for Z^{xi + ri) at 
the independent point x[ — Xi+ri. By substraction, we obtain an equation for the differences 
AZf' — {Zfy — Zf (here and in the following "primed" variables are intended as calculated 
on the point x[). Using the hypothesis of independence of points a;^ and Xi with respect to 
derivatives, namely di{Z^y — d'^Zp — (where 9- represents derivative with respect to x[), 
we get 

dtAZf + Z^d'^AZf = -{dl + d,)AP + {dl' + dl) [A±AZ+ + X^AZr] (2) 
{AP = n' -tt). By adding and substracting the term Z^d'^AZf to © we obtain 

dtAZf + AZ^d'^AZf + Z^{d'^ + a„)AZ± = -{d[ + di)AP + 

+ {dl: + dl)[\^AZt + X^AZr] (3) 

We are seeking for an equation for the second-order correlation tensor {AZfAZ^) related 
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to pseudo-energies. In fact, in a more general approach one should look at a mixed tensor, 
namely (AZ^AZj), taking into account not only both pseudo-energies but also cross-helicity 
{Z+Z-) and {Z-Z+). However, using the DIA closure by Kraichnan, it is possible to show 
that these elements are in general poorly correlated [13j . Since we are interested in the energy 
cascade, we limit ourself to the most interesting equation that describes correlations about 
Alfvenic fluctuations of the same sign. To obtain the equations for pseudo-energies we multiply 
equations ([3|) by AZ^ , then by averaging we get 

dtiAZ^AZf) + {AZ^dU^zfAZf)) + + 9„)(AZ±AZ±)) = 

= - {AZf{dl + d,)AP + AZf{d'^+dj)AP) + 

+ A±(AZ±(9^' + dl)AZ+) + A±(AZ±(a^' + dl)AZ+) + 

+ X'^iAZfidl' + dl)AZ7) + X^{AZf{dl + dl)AZJ) (4) 

If we consider local homogeneity we have 



d{x'^ - ra) dra 

when applied to difference quantities, so that the nonlinear term, using incompressibility, 
becomes 

(AZ^d'JAZfAZf)) = A(AZT(AZ±AZ±)) 

Note that in eq. ([T|) kinematic viscosity are not assumed equal to magnetic diffusivity, 
and this generates a coupling between Z^ and Z^ not only in the nonlinear term but also 
in the dissipative term. We exclude these couplings by making here the usual symplifying 
assumption, that kinematic viscosity is equal to magnetic diffusivity, = = v. Then, 
by using the independence of derivatives with respect to both points and using the local 
homogeneity hypothesis, the dissipative term becomes 

.((a- + dl){AZfAZf)) = 2.^{AZfAZf) - ^^(e,=^r„) 

where we defined the average dissipation tensor 

6± = i.{{d^Zf){d^zf)) (5) 
Using these equations in (|4]), we finally obtain the equation 



dtiAZfAZf) + —{AZ^iAZfAZf)) 



Or, 
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The first and second term on the r.h.s. of tlic last equation represent respectively a tensor 
related to large-scale inhomogeneities kij = {Z^ {d'^ + da){AZf AZ^)) , and the tensor related 
to the pressure term = {AZf{d[ + 9i)AP + AZf (9j + dj)AP). Equation (O is an exact 
equation for anisotropic MHD equations that links the second-order complete tensor to the 
third-order mixed tensor through the average dissipation rate tensor. 

Using incompressibility and independence of derivatives with respect to both points, the 
first term on the r.h.s. can be written as Aij = {d'^ -\- da){Z^ {AZf AZ^)) which vanishes for 
a globally homogeneous situation, because in this case di{) = 0. The pressure term is more 
complicated to be managed. Using independence of derivatives and local homogeneity we get 



(AZfd.AP) = (9,[AZ±AP]-(9,z±)AP) = 
= -{d'[AZ^AP])-{{d,zt)AP) 



(7) 



from which 



n,, = {[d'^iztY - d,zf] AP) + {[d[{zf)' - d^zf] AP) 



(8) 



Then the diagonal terms of the tensor containing the pressure vanish. In fact summing over 
indices eq. ^ yields [d[{zf)' — dizf] which is zero for local homogeneity and incompressibility. 
This means that, assuming global homogeneity and incompressibility, the equation for the 
trace of tensor can be written as 



dt(\AZl 



±|2\ 
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This expression is valid even in the anisotropic case, that is fields depends on the vector r^. 
Moreover by considering only the trace, we ruled out the possibility to investigate anisotropies 
related to different orientations of vectors within the second-order moment. It is worthwhile 
to remark here that only the diagonal elements of the dissipation rate tensor, namely are 
positive defined, while in general the off-diagonal elements can be in principle also negative. 
For a stationary state the equation ^ can be written as the divergenceless condition of a 
quantity involving the third-order correlations and the dissipation rates 



d 
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±|2\ 
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from which we can obtain the Yaglom's relation by projecting equation (jlOp along the longi- 
tudinal — rSr direction. This operation involves the assumption that the flow is locally 
isotropic, that is fields depends locally only on the separation r, so that 



d_ 

dr 



{AZ^\AZ, 



±|2\ 



- 2v 



d_ 

dr 



\AZt\^) + -etr 



= 



(11) 



The only solution that is compatible with the absence of singularity in the limit r is 



{AZ^\AZt\^)=2v^^{\AZ 



(12) 



which reduces to the Yaglom's law for MHD turbulence as obtained by Politano and Pouquet 
[lOj in the inertial range when — > 



(AZ,T|AZ±|2)=--4r 



(13) 
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Finally, in the fluid-like case where zf = = Vi we obtain the usual Yaglom's law (Aw^ | p) = 
—4/3 (er) (e being the usual dissipation rate) which immediately reduces to the Kolmogorov's 
law (Awj?) — — 4/5 (er) in the isotropic case where {Avr^Vy) — {Avr^vD = l/3(Awj?) (as- 
suming the separation r along the streamwise direction x). 

Even if eq. ([T^ remains formally valid only for isotropic MHD, this cannot completely 
solve the problem of the energy cascade in MHD. The tensor liij is zero only when local 
anisotropy is assumed. In fact by calculating the divergence with respect to the index i of 
Iljj, assuming independence of derivatives, we get 

d,{^zf{d,P - d]P') = (AZ±a.a,P> = 9.(AZ±9,P) = -(9,'AZ±9,P) = (14) 

By symmetry the divergence with respect to j also vanishes, that is dfjUij = 0. By using for 
By the isotropic formula for a generic tensor [18] 



n,,(r) = [nn(r) - B„„(r)] ^ +n^^{r)S, 



r 
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it can be easily shown that B^j = by local isotropy [TS] . 

Reduced MHD. - As a different approach, let us consider the Reduced MBD approxima- 
tion (RMBD) [IHllSn] which is valid under the hypothesis of a strong guide magnetic field. 
In this approximation, the dynamics along the parallel and perpendicular directions (with 
respect to the average magnetic field) are disentangled. Then we can distinguish between 
derivatives for perpendicular {d/dx±) and parallel {d/dx\\) coordinates. In terms of Elsasser 
variables, the RMHD approximation reads [H] 

^ + Z^^±^ = -^ + A±^ + A^^ (15) 

dt ^ dx± dx\\ dx±^ dx\ dx\ 

where Z^(x\\,x±_) is the perpendicular component of the Elsasser variables. From this equa- 
tions, by performing the same calculations as before, and by defining the separations along 
the perpendicular r±^ and parallel r\\ directions respectively, we obtain the following Karman- 
Bowarth relation for the stationary state 

(AZJ |AZf = 2^-^(|AZ^ p) - 2e±r^ + - A(|AZf ^drl (16) 

From this last equation it is evident that, in the limit of vanishing viscosity i/ — > 0, in general a 
Yaglom's relation cannot be derived. This can only be the case by assuming that the average 
pseudo-energies are almost constant along the parallel direction. So, if the third term on the 
r.h.s. is zero, we can derive a Yaglom's relation between the mixed third-order correlation 
term as a linear function of the transverse scale rj_ 

(AZj|AZfp)^-24r^ (17) 



Even if this is only an approximate relation, it is worthwhile to note that, in this last case, 
by assuming the scaling AZj^ 

I 2/3 

AZj^ ~ , which imediately le 
anisotropic MBD turbulence [8l. 



by assuming the scaling AZ^ ^ AZ^, relation (fT7|) is compatible with the scaling law 
AZ^ ~ r']^^ , which imediately leads to the Kolmogorov spectrum E{k±) ^ k^/"^ predicted for 
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Conclusions. - To conclude, we reviewed the derivation of a general Yaglom's equation 
for MHD turbulence where two fields are coupled. When this equation is satisfied a turbulent 
cascade is at work. The most general equation ^ that relates the third-order mixed tensor 
to the dissipation rate tensor, is valid in the anisotropic and nonhomogeneous case. By using 
homogeneity, we derive the equation ([T0|) . that is valid in presence of anisotropy, while to 
derive the usual Yaglom's law ([T3|) the local isotropy assumption is required. Moreover, the 
tensor containing the pressure has zero diagonal elements, but is completely zero only when 
local isotropy is assumed. As far as the local isotropy assumption is considered, while for 
usual fluid flows the return to isotropy at small-scales is assured, even if they are anisotropic 
at large-scale [H], in general th MHD flows cannot return completely to isotropy at small 
scales P^UTT] . This in turn means that if we want to fully investigate the turbulent cascades 
in anisotropic MHD flows, the off-diagonal elements of the third-order mixed tensor cannot 
be disregarded. Of course the tensor term related to pressure cannot be eliminated, while 
if we cannot assume local isotropy the relation pH)) cannot reduces to the usual Yaglom's 
law obtained in Ref.s |iOj and investigated experimentally in Ref.s [Ullll]. Of course, since 
the tensor term Hij cannot be calculated using solar wind data, we cannot have any feeling 
of the relative importance of this term and the term containing the dissipation rate tensor 
in anisotropic MHD turbulence. High-resolution numerical simulations for anisotropic MHD 
turbulence should be used as a first approach. 

* * * 

We acknowledge useful discussions with W.H. Matthaeus. 
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